Thermal activation of vortex- ant i vortex pairs in quasi-2D Bose-Einstein condensates 
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Here we show, by performing ab initio classical field simulations that two distinct superfluid 
phases, separated by thermal vortex-antivortex pair creation, exist in experimentally producible 
quasi-2D Bose gas. These results resolve the debate on the nature of the low temperature phase(s) 
of a trapped interacting 2D Bose gas. 



The phenomena of superconductivity and superflu- 
idity are striking manifestations of the role played by 
quantum statistics at low temperatures. Altering the 
temperature or effective dimensionality may radically 
change the physical properties of quantum degenerate 
systems. A well known consequence of this is that, un- 
like in 3D, there is no Bose-Einstein condensation (BEC) 
for a homogeneous 2D ideal-gas in the thermodynamic 
limit at any finite temperature Nevertheless, the 

Berezinskii-Kosterlitz-Thouless (BKT) vortex binding- 
unbinding phase transition allows the emergence of su- 
perfluidity in 2D systems [^^|^. This superfluid tran- 
sition has been experimentally observed in liquid helium 
thin films 0, superconducting Josephson-junction ar- 
rays 0, and in spin-polarized atomic hydrogen 0|. Al- 
though weak particle interactions alone are not sufficient 
to change the situation, an external confinement modi- 
fies the density of states in such manner that the criti- 
cal point of BEC is elevated to finite temperatures 0- 
Therefore it is not certain a priori whether the transfor- 
mation from normal to superfluid in such systems is BEC 
or BKT-type transition. 

Significant advances in trapping ultra-cold atoms, and 
in particular the use of optical lattices, has made pro- 
duction of quasi- 2D quantum gases experimentally fea- 
sible [Tol Ull IT^ IT3I ll^ . Low temperature ordering in 
these systems leads to the competition between BEC and 
BKT-type transitions. Long wave-length phase fluctua- 
tions have been predicted to destroy the phase coherence 
of a pure condensate at low temperatures leaving the 
system in a state of so-called quasi-condensate, however, 
the structural details of such a state have remained under 
speculation [UlllllllillllliSllllll. In this Let- 
ter we provide quantitative evidence for the existence of 
a vortex-free condensate phase and a BKT-type phase 
manifested by the presence of vortex-antivortex pairs, 
using a formalism that allows direct observation of the 
transition between these two distinct superfluid phases. 
Furthermore, we show that the phase defects observed in 
the experiment by Stock et al. [lO| are readily explained 
in terms of spontaneous activation of vortex-antivortex 
pairs in the temperature regime of the predicted BKT- 
type phase. 

The statistical probability, p, for the excitation of a 
vortex-antivortex pair in a phase coherent Bose-Einstein 
condensate is proportional to the Boltzmann factor, 
^-F/kuT ^ where F is the free energy cost of such topo- 
logical excitation, fee is the Boltzmann constant and T 



is the temperature. The condition p = 1 may be used 
to estimate the critical temperature, Tc, separating the 
vortex-free and vortex-pair phases. Considering a single 
pair in the vicinity of the 2D harmonic trap centre, the 
hydrodynamic approximation yields [l5l | 
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where, /3 = 27r'^ft^(ln 2-|~e)^/3TO, A^o is the number of par- 
ticles of mass m in the condensate phase, g is the strength 
of the particle interactions, and Tg is the critical temper- 
ature for an ideal-gas to undergo Bose-Einstein conden- 
sation. The constant e is the core energy per particle of 
a vortex pair, expressed in the trap units. 

In order to numerically study the low temperature 
phase diagram of a trapped quasi-2D Bose gas and the 
validity of Eq. we divide the system into classical 
and so-called incoherent regions depending on the mode 
occupation of the single particle eigenstates. The classi- 
cal region consists of the highly occupied modes that are 
described by the projected Gross-Pitaevskii equation 
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where Vc^t denotes an external potential and the projec- 
tor V serves to restrict evolution of the classical field, 
^'(r,f), to within its subspace [2^|3. The incoherent 
region consists of the modes of low occupation which are 
treated using the semi-classical Hartree-Fock approxima- 
tion. The classical and incoherent regions are taken to be 
in equilibrium with each other. The former is most con- 
veniently represented using the harmonic oscillator eigen- 
states with energies below certain cutoff energy, Ed- This 
energy cutoff needs to be chosen such that the mean oc- 
cupation of all the modes in the classical region is > 1. 
The normalization of '^(r,t) determines the number of 
particles Nc\ > Nn in the classical region. To simulate 
the experiment |10| , we concentrate on a single quasi-2D 
site of their optical lattice and perform simulations em- 
ploying a 3D harmonic trap with respective radial and 
axial frequencies {uj±,uJz} = 2tt x {106,4000} Hz. This 
choice corresponds to the shorter lattice period in the 
experiment. The time and spatial scales in our calcula- 
tions are Iq ~ 27r/wj_ and ao = ^Jh/2muj_]_, respectively. 
Depending on the total particle number, A^, the energy 
cutoff is in the range, E^i = (32 — 56) fiw_L, resulting in 
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FIG. 1: Low temperature phase diagram for quasi-2D Bose 
gas. The markers denote diflerent vortex-pair excitation prob- 
abihties p = (o), < p < 0.25 (®) and p > 0.25 (•) inside 
a circular region whose diameter equals the Thomas-Fermi 
radius. The lines are analytical estimates using Eq. Q with 
e = (— ) and e = 1/2 (--). 

the non-pcrturbativc treatment of several hundred lowest 
modes of the many-body system. 

Each simulation proceeds as follows. The initial mode 
distribution of the classical field, ^'(r, i), is first allowed 
to equilibrate for 10 trap periods, i.e. 10 to, after which 
A^o, N and the fractional temperature, T/Tq, are ex- 
tracted by averaging over 500 microstates, sampled uni- 
formly over another 10 to- The reference ideal-gas con- 
densation temperature, Tq, is obtained by exact summa- 
tion of the partition function for the quasi-2D harmonic 
oscillator. Subsequently we numerically detect vortex 
pairs by locating the phase singularities in the classical 
field, ^'(r, t). To measure the vortex-pair excitation prob- 
ability, p, we only include vortices observed inside a circle 
whose diameter is equal to the Thomas-Fermi radius for 
a pure condensate of Nq particles. Thus we evaluate, p, 
by counting the average number of vortex pairs in the 
500 microstates sampled for each simulation. 

The computed superfluid phase diagram for a trapped 
quasi-2D Bose gas is shown in Fig. ^ The insets show 
typical planar density distributions in the two superfluid 
phases. We have categorized the simulated systems ac- 
cording to the measured vortex-pair excitation probabil- 
ities p = (o), < p < 0.25 (®) and p > 0.25 (•). The 
lines are plotted using Eq. ^ with (--) and without (— ) 
the core contribution e. Two distinct phases — vortex-free 
and vortex-filled — are evident in the diagram. The an- 
alytic result provides a fair estimate of the critical tem- 
perature for the transition between the two superfluid 
phases. The temperature range in which phase defects 



FIG. 2: Internal energy for = 1 x 10* particles (main frame) 
and the vortex excitation probability for a range of A'^ (inset) 
as functions of the fractional temperature. The solid line is 
a fourth order polynomial fit to the data. The dashed line is 
the ideal-gas result. 

were observed in the experiment is consistent with that 
of the BKT-type phase in Fig. ^ The slight upward 
tendency in the computed transition temperature with 
increasing N is probably due to axial degrees of freedom 
becoming thermally excited, causing the system charac- 
teristics to approach three dimensionality. 

Figure [5| shows the internal energy per particle of 
the function of the fractional temperature for 

N = 10^ particles. The computed vortex-pair excitation 
probability, p, corresponding to the data in Fig. ^ is dis- 
played in the inset. The solid line in the main figure is 
a fourth order polynomial fit to the data. The interac- 
tions impose an overall increase in the internal energy 
in comparison to the ideal-gas result (--). This energy 
characteristic alone is too coarse to reveal latent heat at 
the vortex-pair creation transition in this inhomogeneous 
mesoscopic system. However, as is evident in the inset, 
the measure of vorticity provides a useful order parame- 
ter to characterize the transition. 

At the lowest temperatures no vortices are present in 
the system. They first emerge in the low density regions 
of the cloud and are only able to nucleate closer to the 
trap centre as the temperature increases. We note that 
at intermediate temperatures the phase coherent central 
condensate and outer shell of phase fiuctuating superfluid 
co-exist, while at even higher temperatures the whole 
cloud turns critical with vortex pairs able to form and 
annihilate everywhere. Ultimately, as temperature in- 
creases, superfluidity is lost when the system crosses the 
BKT-type transition to the normal state. Thus the exter- 
nal confining potential, which fundamentally allows BEC 
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to exist in 2D. also leads to a spatial dependence of the 
critical temperature for the spontaneous vortex-pair cre- 
ation. The lifetime of vortex-pairs is dependent on the 
local superfluid density, being shorter in the high density 
central regions. 

In the experiment a set of quasi-2D systems were re- 
alized using an optical lattice, and were subsequently al- 
lowed to expand in free space and interfere with each 
other. The presence of phase defects were inferred from 
the distortion of the interference patterns. In order to 
fully explain those observations, we have simulated this 
interference procedure. We construct a 3D wavefunction 
corresponding to a lattice of quasi-2D systems, where 
the individual components are taken as independent mi- 
crostates from a single classical field simulation. The 
lattice parameters are chosen to match those used in the 
experiment. Subsequently, we propagate thus prepared 
wavefunction according to a trap-free time-dependent 
Gross-Pitaevskii equation. Upon release from the exter- 
nal potential, mean-field effects on the evolution rapidly 
decay and the momentum distribution freezes, allowing 
us to infer the asymptotic spatial interference patterns. 
In practise they are obtained by Fourier transforming the 
spatial wavefunctions after a 0.5 io time-of-flight. 

Figure |2Ja) shows a result of such numerical interfer- 
ence experiment for the two classical fields whose pla- 
nar density distributions are shown in Figs.l^Jb) and (c). 
The characteristic zipper pattern in the interference is 
predominantly caused by the phase singularity associ- 
ated with the isolated vortex in the centre of Fig. |2Ib). 
Despite the transient imbalance in the number of oppo- 
sitely charged vortices in Fig. I^fb), the value of angular 
momentum remains zero by virtue of the spatial vortex 
distribution. At low temperatures only straight fringes 
are produced in the interference patterns, manifesting the 
global phase coherence of the condensates. In the pres- 
ence of vortex pairs, zippers and more subtle structures 
emerge revealing the underlying phase fiuctuations. 

We have verified that the zipper pattern observed in 
the experiments is indeed produced by an isolated vor- 
tex and that closely bound vortex-antivortex pairs do 
not significantly alter the interference pattern. However, 
such isolated vortices are remnants of vortex pairs which 
have dissociated, since conservation of angular momen- 
tum should prevent formation of unpaired vortices. In 
the vicinity of the vortex-pair creation transition only 
tightly bound vortex pairs exist but at higher temper- 
atures, deeper in the BKT-type phase, seemingly inde- 
pendent vortices are abundant in the system since the 
energy cost of unbinding an existing pair is comparable 
to the cost of creating a new pair. We have verified such 
dynamic behaviour by visually observing the creation of 
vortex pairs and their subsequent dissociation and/or an- 
nihilation in simulations. The isolated vortices are most 
clearly identified by the interference pattern they pro- 
duce, whereas the lack of resolution in this type of detec- 
tion method obscures the observation of tightly bound 
vortex pairs. This, together with the relatively narrow 



temperature range of the BKT-type phase for these sys- 
tem parameters, could explain the rather low probability, 
~ 10%, for observing a clear signature of phase defects in 
the experiments. Therefore, further development of de- 
tection methods may be required to observe the onset of 
the vortex-pair creation (i.e. BEC-BKT) transition. The 
probability for detecting phase defects could perhaps be 
enhanced by imparting angular momentum to the system 
to artificially separate the paired vortices. 
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FIG. 3: Interference pattern (a) produced by two independent 
classical fields (b) and (c) at temperature T = 0.86 To. The 
relevant particle numbers are A^ci = 3.0 X 10^ and Af = 4.0 x 
10'*. The zipper structure in (a) is the telltale signature of 
the phase singularity associated with the central vortex in 
(b). The locations of vortices and antivortices are marked by 
4- and — signs, respectively. 

Vortex-antivortex pairs are the only phase defects ob- 
served in the simulations. Other possibilities, such as 
dark solitons, can be excluded from the picture since 
their cost of energy, F , is much greater than that of a 
vortex(pair). And even if excited, oppositely travelling 
solitons would rapidly decay into vortices via the snake 
instability. However, at very high temperatures the vor- 
tices are occasionally seen to arrange themselves into a 
chain and percolate throughout the cloud in anticipation 
of the transition to a normal state. 

To conclude, we have performed ah initio classical field 
simulations of quasi-2D Bose-fields and have character- 
ized the low temperature phases for such systems over 
a wide parameter range. We have provided strong evi- 
dence supporting the view that the BKT-type phase was 
observed in the recent experiment by Stock et al. |lOl |. 
Moreover, we have demonstrated the limited sensitivity 
of the experimental detection method, based on pairwise 
interference of individual clouds, to tightly bound vortex- 
antivortex pairs. Our results further vindicate the appli- 
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cability of classical field methods to the study of ultra- 
cold atomic gases. An additional benefit of our approach 
is that it allows for a direct study of the proliferation of 
vortices in the normal to superfluid (BKT-type) phase 
transition. 
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